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ABSTRACT 


ill 


Since  the  pioneering  work  of  Giaever  (i960  a,b), 

the  electron  tunneling  technique  has  been  successfully  applied 

■ 

to  the  study  of  energy  gap,  quasi -particle  density  of  states 
and  phonon  effects  in  superconductors. 

Scalapino  et  al  (1965)  predicted  a  peak  in  the 
effective  tunneling  density  of  states  in  lead  at  non-zero 
temperatures  due  to  recombination  effects.  Using  the  electron 
tunneling  technique,  we  have  obtained  results  for  the  varia¬ 
tion  with  voltage  of  the  normalized  dynamic  conductance  of 
aluminum/aluminum  oxide/lead  tunnel  junctions  which  agree  well 
with  their  published  work  and  give  definite  evidence  for  the 
occurrence  of  the  recombination  phenomena.  The  energy  gap 
for  lead  was  found  to  be  2A(0)p^  =  (2.68  ±  0.06)  mev  which 
is  also  in  agreement  with  published  results. 

During  the  course  of  the  present  work,  an  interpola¬ 
tion  formula  for  calibration  of  a  germanium  thermometer  using 
only  three  known  temperatures  was  obtained  in  the  range  of 
temperatures  between  7*2°  K  and  0.52°  K.  The  interpolated 
temperatures  are  accurate  to  ±  0.05°  K  above  2.5°  K  and  to 
nearly  1%  below  2.5°  K. 
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CHAPTER  I 


INTRODUCTION 


A .  Motivation 

This  thesis  reports  an  extension  of  work  on  the 
electron  tunneling  in  superconductors  which  was  initiated  in 
this  laboratory  by  Adler  (1963). 

According  to  the  Bardeen-Cooper-Schrief fer  theory 
of  superconductivity  (henceforth  called  the  BCS  theory)  a 
forbidden  energy  gap  appears  about  the  Fermi  energy  in  the 
density  of  electron  states  when  a  metal  turns  superconducting. 
The  electron  tunneling  technique  by  which  electrons  of  varying 
energies  may  be  injected  into  a  superconductor  provides  a 
direct  and  simple  method  of  observing  this  energy  gap.  One 
can  obtain  not  only  the  magnitude  of  the  energy  gap  but  also 
the  variation  of  the  density  of  states  with  energy  outside  the 
gap  in  a  superconductor.  Much  smaller  variations  in  the 
density  of  states  above  the  energy  gap  have  also  been  studied 
using  the  electron  tunneling  technique  and  have  been  identi¬ 
fied  as  due  to  the  predominant  phonon  energies  of  the  metal. 
Thus,  while  providing  a  testing  ground  for  the  BCS  theory  and 
its  modifications,,  this  simple  technique  has  also  led  to 
investigations  of  the  phonon  spectrum  of  several  superconductors 
and  holds  promise  of  providing  a  great  deal  of  information  on 


1 


•  ;  ■  :S  tal  **"  8-xe  k.  ifto  ■.  k  ii  Snil  .nnu  to itosie 


t  at  .'Moot -agaa  urutf  lactam  s  harfw  awte3a  noxcfoaia  lo  ^xanaft 
a',iV  V  UC  *****•!•  rtoitiw  *  a.  pinrfoa.1  s^ixacawi  no^oala  artT 
•'•  £  ?J  ^Woulmoowqua  a  o^ni  besostnl  »d  ^  eaXyxena 

**  Vi  09  rWJw  s^Bia  *  m  en.t  lo  nott*lMv  -,r  ■ 

ri  ,tJ  8  oti£*  9V£"  **8  Wisne  ,rii  svodB  ao3Ma  lo 

nsbi  ne  d  bns  -o.  snIXaanwd  ncudoaie  art3  gniau 

3  x  3  «*  arid  ,ol  bnuoxg  a  gr,  ibxvoxq  aXJ:<*r  ,?:Jrt, 

•“  0  rr onorlq  erfcf  lo  anoldBSiiaevnI 

' 


2 


the  properties  of  superconductors  in  general  e.g.  the  double 
particle  tunneling  etc.  (of.  chapter  III). 

The  specimens  used  in  the  present  work  consisted 

of  a  thin  insulating  layer  sandwiched  between  two  thin 

metallic  films.  These  specimens  are  ideally  suited  to  low 

temperature  work  because  of  their  small  size.  The  current 

through  the  insulating  layer  was  measured  as  a  function  of 

the  voltage  across  it.  If  the  potential  difference  across 

the  layer  is  small.,  the  current  that  tunnels  through  it  is 

directly  proportional  to  the  potential  difference  as  long  as 

the  density  of  states  in  the  two  metals  is  constant  over  the 

applied  voltage  range  (cf.  Fisher  and  Giaever  1961).  When 

one  or  both  of  the  metallic  films  turn  superconducting,,  the 

density  of  states  will  vary  rapidly  in  a  small  voltage  range 

and  the  tunneling  current  will  no  longer  vary  linearly  with 

the  applied  voltage.  This  non-linearity  is  related  to 

variation  in  the  density  of  states  in  the  superconductor.  In 

a  metal/insulator/superconductor  system  (henceforth  referred 

to  as  a  N-S  system)  the  normalised  dynamic  conductance  cr, 

di 

which  is  the  ratio  of  the  dynamic  conductance  (^)  when 

ns 

one  of  the  metallic  members  of  the  tunnel  junction  is  super- 

conducting  to  the  dynamic  conductance  (tv-)  when  both 

av  nn 

members  are  normal,  is  measured  as  a  function  of  applied 
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voltage  v.  A  graph  showing  a  as  a  function  of  v  directly 
displays  the  variation  in  the  normalised  density  of  states 
with  energy  in  the  superconductor.  The  normalised  density  of 
states  in  a  superconductor  is  the  ratio  of  density  of  states 
when  a  metal  is  superconducting,  to  the  density  of  states 
when  it  is  normal. 

The  purpose  of  the  present  work  was  to  look  for  a 
peak  in  the  effective  tunneling  density  of  states  in  lead 
predicted  by  Scalapino  et  al  ( 1 965 )  to  be  present  at  non-zero 
temperatures  because  of  recombination  processes  that  had 
been  neglected  in  previous  theoretical  work. 
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CHAPTER  II 

MICROSCOPIC  THEORIES  OF  SUPERCONDUCTIVITY 
A .  The  BCS  Theory 

Superconductivity*  first  observed  by  Kamerlingh  Onnes 
in  1911*  is  one  of  the  most  fascinating  properties  of  metals. 

It  was  only  in  1957  that  the  first  successful  microscopic 
theory*  which  explains  most  of  general  features  of  supercon¬ 
ductivity*  was  proposed  by  Bardeen*  Cooper  and  Schrieffer 
(1957  a, b) . 

The  basic  difficulty  in  constructing  a  theory  of 
superconductivity  is  that  the  difference  in  energy  between 
normal  and  superconducting  phases  --  the  condensation  energy  -- 
is  extremely  small.  Hence  one  has  to  pick  out  the  terms  in 
energy  which  differentiate  between  two  phases  and  calculate 
this  difference  directly. 

BCS  theory  has  its  origin  in  work  of  Cooper  (1956) 
who  was  able  to  show  that  if  there  is  a  net  attractive  force 
between  two  electrons  excited  slightly  above  the  Fermi  surface* 
they  can  form  a  real  bound  pair  ha  vi^g-^^Q-cadu^^d:rwa^e~-f ane-t ion- . 
The  energy  of  a  system  containing  bound  pairs  is  less  than  the 
ground  state  of  the  Fermi  sea.  Thus  the  ground  state  of  the 
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Fermi  sea  is  unstable  against  the  formation  of  such  Cooper 
pair  states  and  the  total  energy  would  be  reduced  from  that  of 
the  "normal”  state.  The  BCS  theory,  then,  provides  a  formalism 
for  dealing  with  a  state  in  which  a  large  fraction  of  the 
electrons  exist  in  such  bound  pairs. 


The  simplest  way  to  describe  the  attractive  force 


between  a  pair  of  electrons  is  based  on  virtual  phonon  exchange 
between  them.  The  matrix  element  for  scattering  electrons  from 
states  Tc  and  IcT  to  Ic  +  J?  and  IcT  -  It  by  exchange  of  a  phonon  is 

2 

(2.1) 


214  1  "g  g+  It 

" 


where  Mj*-  ^  +  ^  is  the  matrix  element  of  the  electron-phonon 

interaction,  di  is  the  energy  of  the  phonon  involved  and 

Q. 

€£•  is  the  Bloch  state  energy  measured  from  the  Fermi  level. 
The  interaction  is  attractive  (negative)  if 


elt  sIc+lf 


<  *«  “5  *  kB  % 


(2.2) 


where  is  the  Boltzmann  constant  and  0^  is  the  Debye 

temperature.  In  this  range  of  energies,  the  magnitude  of 

expression  (2.1)  is  roughly  equal  to 

2 
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To  this  we  have  to  add  the  repulsive  Coulomb 
.attraG-ti-offr  for  which  the  matrix  element  would  be 


2 

4  7T  e 


+  A 


2 


(2.4) 


where  A  is  a  screening  constant  and  e  denotes  the  electronic 
charge . 

The  criterion  for  formation  of  bound  pairs  and 
superconductivity  is  that  the  attractive  interaction  dominates 
the  screened  Coulomb  interaction.,  This  criterion  can  be 
expressed  in  the  form 


-  V  = 


<- 
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1  Mg,  g  +  g 1 

q 


4  it  e^  \ 

2  / 

K  +  A  x 


(2.5) 


where  the  average  has  been  taken  over  electron,  states  near 
the  Fermi  surface  which  satisfy  Eqn.  (2.2). 


Since  electrons  follow  Fermi -Dirac  statistics,  the 
matrix  elements  for  all  possible  interactions  which  transfer  two 
electrons  from  any  two  Ic  values  to  any  two  others  satisfying 
Eqn.  (2.2)  alternate  in  sign  depending  on  the  occupation  number 
of  other  states  which  remain  unchanged  in  the  transition: 
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<  k^,k^  .  .  .  Jc+It,  IcT  -It  |  V  |  'Icp,  .  .  .  Ic,  Tc!  > 

=  ±  <  Ic+It,  IcT  -It  |  V  |  lx,  IcT  > 

The  step  taken  by  BCS  is  to  associate  all  the  Ic 
values  in  pairs  (k,!?’ )  and  require  that  either  both  or  neither 
member  of  a  pair  be  occupied.  This  will  result  in  a  single 
sign  of  the  matrix  element. 

The  momentum  must  be  conserved  in  scattering  pro¬ 
cesses.  The  maximum  number  of  possible  transitions ,  hence 
the  lowest  energy,  can  be  obtained  only  when  the  momentum  of 
each  pair  is  the  same  so  that  any  pair  can  be  scattered  into 
any  other  pair.  Thus  (lx  -h  Tc T  )  =  "c  should  be  same  for  all 
the  pairs.  For  the  ground  state,,  the  phase  space  available 
for  transitions  becomes  maximum  when  the  total  momentum 
(T  =  0.  Further,  as  the  exchange  terms  tend  to  reduce  the 
(negative)  interaction  energy  for  pairs  of  parallel  spin,  it 
is  most  probable  that  the  paired  state  should  contain  opposite 
spins.  Thus  the  superconducting  ground  state  can  be  expressed 
entirely  in  terms  of  paired  states  (let  ,  -IcJ')  which  are 
occupied  simultaneously.  The  condensation  energy  of  super¬ 
conductivity  will  arise  from  interaction  terms  which  scatter 
a  given  pair  of  electrons  from  one  such  paired  state  to  another 
similarly  paired  state.  Thus  basically  a  two  body  correlation. 
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with  strong  preference  for  zero  momentum  singlet  pairs,  is  the 
cause  of  superconductivity „ 


The  Hamiltonian  corresponding  to  this  condensation 
energy  -  often  referred  to  as  the  reduced  Hamiltonian  -  of  the 
system  will  be 


H 


red 


=  T  +  V 


red 


=  22 
Ic 


'k 


* 

by*  by* 

k  k 


2 


V- 


kkt 


* 

by*.  by* 

kT  k 


(2.6) 


Here  T  is  the  kinetic  energy  term,  is  the  interaction 

potential  acting  to  scatter  a  Cooper  pair  from  states  (Tct,  -Ic±) 
to  (Tc T T  ,  -  It’i)  by  an  interaction  strength  ,  b^*  and 

b^*  are  pair  creation  and  annihilation  operators  defined  by 


and 


*  *  * 
tc  =  °icT  C-Hlr 

bIf  =  °~faCZT 


(2.7) 


Cr>y  and  Ct*  /  are  the  creation  and  annihilation  operators 
k/  O'  k/g 

respectively  for  electrons  in  Bloch  states  with  crystal  momentum 
£  and  spin  index  a„  They  satisfy  the  Fermi  anti commutation 
rules  and 

■  =  ct /a  ct/d  (2‘8) 


is  the  number  operator  whose  expectation  value  gives  the  proba¬ 
bility  that  state  k/cr  is  occupied. 
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BCS  assume  that  the  probability  of  a  specific 
configuration  of  pairs  is  given  by  the  product  of  probabilities 
of  each  individual  pair  state  being  occupied.  Thus  the  ground 
state  function  is 


T 

o 


k 


Ut-*  + 

■  k 


*  V 

Vr^bv*  ) 

k  k  ' 


0 

o 


l « 


i/i  V*  SA.IAYVV 

2  ^ 
where /  v  ^  is 


and  u  ^  is  the 


the  probability  that  the  k 
probability  that  it  is  not 


th 


pair  is 


occupied. 


(2.9) 

occupied 

Clearly 


2 

U  g  + 


2 

V  1? 


=  1 


(2.10) 


The  energy  corresponding  to  the  reduced  Hamiltonian 

H  ,  for  the  state  (2.9)  is 
red  ' 


WC  =  2  l  ef  v  T?  -  TV,  VB?'  u£  v£  U^’  V^’  • 

°  k  kk T 


(2.11) 


Treating  v-£  as  a  variational  parameter,  the  value 
of  v^  which  gives  minimum  value  of  Wc  is 


and 


v 


u 


=  1  (i  -  ^  ) 


k 


K 


k 


2 


—  —  ( i  +  — 

2  1  E 


E. 


k 


(2.12) 


k 


Tc 


where 


E 


k 


=/  +  A2rir 

v  k  k 


(2.13) 
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It  follows  that 

=  2  AU' 

k  ^  - 

k'  2  E**, 

kT 


(2.15) 


BCS  make  the  following  assumptions: 

(a).  (2.1)  is  taken  to  be  constant  given  by  (2.3)  in  the  range 

where  |  €-*  [  and  I  e]£  +  it  I  are  less  than  'fi  uw 
and  zero  elsewhere. 


(b) .  An  average  value  of  (2.4)  is  used  and  is  taken  to  be 
constant  in  the  same  range  of  energies  as  in  (a)  and 
zero  elsewhere. 


~  constant 
k 

otherwise . 


A 


for 


v<  and  zero 

s  q 


(d)  .  The  density  of  states  N(e^>)  is  assumed  to  be  constant 
N(0)  in  the  range  |s^  |  ^  'fi 


With  these  assumptions  Eqn.  (2.15)  reduces  to 


-T5 

A _ Q 

"  sinh  (1/  N(O)V) 


(2.16) 


and  the  energy  difference  between  normal  and  superconducting 


states.,  i.e.  the  condensation  energy.,  is 

2 


O 


w  -  w 

s  n 


W  = 
c 


N(0)  (fi  au) 


exp  [2/N(0)v]  -  1 


(2.17) 
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Following  Eqns .  (2.12)  and  (2.13),  the  pair 

occupation  function  in  the  BCS  ground  state  shown  in  Fig.  2.1(a) 

is  seen  to  be  like  a  Fermi  function  with  kJT  ~  .  There 

B  k 

are  few  holes  deep  in  the  Fermi  sea  and  few  highly  excited  pairs. 

In  the  weak  coupling  limit,  i.e.  when  N(0)v«  1, 

the  energy  gap  at  T  =  0°K  is  found  to  be 


2A(0)  =  3c 52  kBTc 

where  T  is  the  transition  temperature.  For  0  <  T  <  T 
the. gap  can  be  numerically  calculated.  The  temperature  depen¬ 
dence  of  the  gap  is  shown  in  Fig.  2.2.  The  excitations  involved 
in  this  calculation  are  discussed  below. 

B.  Elementary  Excitations  and  Density  of  States 

;  ,  theory  of  superconductivity  exactly  equivalent  to 
the  BCS  theory  was  formulated  by  Bogoliubov  (1958)  who  intro¬ 
duced  a  new  set  of  quasi -particle  operators: 


7Ic  “  ul?  CI?  “  VT?  c-t 

and  (2.18) 

*  * 


=  u 


i?  c-i? +  Vir 


This  transformation  preserves  the  Fermi  anticommu¬ 
tation  rules.  The  BCS  ground  state  is  vacuum  for  these  parti¬ 
cles.  These  operators  can  create  and  destroy  excitations  of 
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Figure  2 . 1 


(a)  Pair  occupation  function  in  the  BCS  ground  state. 


(b)  Quasi -particle  energy  in  the  superconductor. 


(c)  Quasi -particle  energy  in  the  normal  metal. 
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Figure  2.2 


Temperature  dependence  of  the  energy  gap. 
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vacuum  in  the  same  manner  as  c^.  and  Cj*  can  respectively 
"create"  and  "destroy"  electrons.  7^  (or  7  always  creates 
an  excitation  with  net  momentum  It  (or  -1?)  and  spin  up  (or  down)  . 
Using  Eqns .  (2.12)  and  (2.13)*  it  is  evident  that  these  excita¬ 
tions  of  the  superconducting  state  are  peculiar  quasi -particles 
which  change  from  being  like  electrons  (e^>  »  A  )  to  being 

like  holes  (e^*  negative  and  |  |  »  A)  as  they  cross  the 

Fermi  surface  at  which  they  are  an  equal  mixture  of  electrons 
and  holes. 


x 

If  f  is  the  BCS  ground  state,  then  7^  ¥  and 
o  °  *  #k  o 

X  X 

7  It  ^  correspond  to  "single"  and  "pair"  excitations 

with  excitation  energies  Ej*  and  2  E^-  corresponding  to* 
respectively*  one  or  both  members  of  a  pair  being  excited. 


Since  excited  single  particles  can  only  be  produced  in  pairs* 
the  minimum  energy  of  an  allowed  excitation  is  2  A.  The 
energy  E^-  of  a  quasi -particle  in  the  superconducting  and  normal 
states  is  shown  in  Figs.  2.1(b)  and  (c) .  These  excitations 
behave  like  independent  fermions. 


Assuming  an  isotropic  energy  gap*  the  density  of 
excited  states  N(E)  in  a  superconductor  is* 
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N(E)  =  N(e)  (dE/de)-1 

where  N(e)  is  the  density  of  states  in  the  normal  metal.  As 

E2  =  s2  +  A2  (2.13) 
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we  have  N(E)  =  N(e) 

-  N(0) 


E 


o 

a"  + 


A 


dA 

de 


E 


using  the  assumptions  (c)  and  (d)  stated  earlier .  As  E 
cannot  have  values  smaller  than  A, 


N(E)  =0  if  |  E  |  <  A 


=  N(0) 


(2.19) 

>  A  . 


Thus  N(E)  is  singular  at  the  gap  edges.  At  energies  well 
above  the  gap,  this  expression  reduces  to  N(0). 


C.  Modifications  to  BCS  Theory 

BCS  theory*  using  only  one  adjustable  parameter 
N(0)V*  compares  surprisingly  well  with  experimental  results. 
However*  anisotropy  in  a  given  metal*  differences  between 
metals  and  structure  in  the  superconducting  density  of  states 
indicate  need  for  improvement  of  this  theory.  In  the  isotropic 
model*  this  is  achieved  by  replacing  V^>t  by  a  kernel  K(e*eT) 
which  has  more  accurate  convergence  properties  at  larger  values 
of  these  variables.  This  leads  to  an  energy  gap  function  that 
depends  on  the  energy  in  a  complicated  manner  and  which 
manifests  itself  in  predictions  of  measurable  quantities. 
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Deviations  from  BCS  theory  are  more  pronounced  in  cases 
where  the  weak  coupling  limit  is  not  applicable. 

In  the  BCS  theory,  the  energy  levels  are  assumed 
to  have  zero  width  so  that  it  is  not  applicable  when  life  time 
effects  are  to  be  considered.  Eliashberg  (i960)  has  included 
life  time  effects  of  these  states  and  phonon  distribution  in 
extending  the  Bogoliubov  theory  with  the  use  of  Green Ts 
function  methods.  Following  this  treatment,  the  energy  gap 
turns  out  to  be  a  complex  function  of  energy  and  has  maxima 
and  minima  at  energies  close  to  0^,  2  ©^  etc.,  (henceforth 
we  will  take  -fi  =  k^  =  e  =  1  where  these  symbols  have  their 
usual  meanings).  Schrieffer  et  al  (1963)  have  shown  that 
Eliashberg  formalism  modifies  the  effective  tunneling  density 
of  states  in  a  superconductor  which  is  now  given  by 


(2.20) 


where  A(E)  is  a  complex  function  of  the  energy. 
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CHAPTER  III 

THEORIES  OF  TUNNELING  IN  SUPERCONDUCTORS 

In  his  pioneering  work,,  Giaever  (i960  a,  b)  studied 
the  tunneling  of  electrons  between  a  superconducting  film  and 
a  normal  one  separated  by  an  insulating  layer  which  was  a  few 
tens  of  angstroms  thick.  Observations  of  current  through  a 
tunnel  junction  as  a  function  of  voltage  with  one  or  both  of 
the  films  in  the  superconducting  state  provide  one  of  the  most 
direct  measurements  of  the  superconducting  energy  gap  and  give 
detailed  information  about  the  density  of  states  and  some  other 
properties  of  superconductors. 

A.  Simple  Model  of  Tunneling 

It  is  well  known  that  a  current  can  pass  between 
two  metals  separated  by  a  thin  insulating  layer  due  to  quantum 
mechanical  tunneling  of  the  electrons,  (cf.  Giaever  and 
Fisher  1961).  The  treatment  has  been  extended  by  Nicol  et  al 
(i960)  and  Giaever  and  Megerle  (1961)  for  the  case  when  one  or 
both  of  the  metals  are  superconducting.  A  good  account  of 
these  theories  has  been  given  by  Adler  (1963). 
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In  this  model,  the  insulating  layer  in  the  tunnel 
junction  is  considered  to  be  a  potential  barrier  for  the 
electrons.  The  transmission  coefficient  for  an  electron  through 

JL. 

this  layer  varies  as  exp  (-th2)  where  t  is  the  insulating 
layer  thickness  and  h  heigh ty.  These  factors *(t  and  h) 

remain  unaffected  by  the  application  of  small  voltages  to  the 
tunnel  junction. 


The  probability  of  transition  of  an  electron  from 
an  unoccupied  state  on  the  left  of  the  barrier  to  a  state 
Tc^,  on  the  right  side  (henceforth  called  metal  1  and  metal  2 
respectively)  per  unit  time  is 


^  I  M  I  2  f±  (1  -  ff)  Nf  (3.1) 

where  M  is  the  matrix  element  for  the  transition,  f^  and 
f  are  the  probabilities  that  initial  and  final  states  are 
occupied,  and  Nf  is  the  density  of  final  states. 


The  net  tunneling  current  across  the  tunnel  junction 
on  application  of  voltage  v  is 


i 


2 


/  2tt  |  M  |2  ^NX(E)  f(E)  (l-f(E-tv))  N2(E+v) 
-00 

-  Ng (E+v)  f  (E+v)  (l-f(E))  N1(E)J  cffil 


(3.2) 


These  parameters  are  reduced  quantities  so  that  the  exponential 
factor  is  dimensionless. 
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where  E  is  the  energy  measured  from  the  Fermi  level  and  f  is 

the  Fermi  function.  M  is  considered  to  be  independent  of  the 

electron  energy  in  the  region  of  interest,  that  is,  when 

\ 

v  «  Ep  being  the  Fermi  energy.  The  density  of  Bloch 
states  near  the  Fermi  level  is  assumed  to  be  constant.  In 
general ^  we  can  write 


NifE)  = 


N.  (0) 
in v  ' 


n.  .(E) 
ij 


(3.3) 


where  i  refers  to  metal  1  or  metal  2  and  J  to  the  normal 
(ri)  or  superconducting  (s)  case  and  N^n(0  )  represents  the 
density  of  states  at  the  Fermi  level.  We  have 


and 


n.  (E) 
in '  ’ 


1 


niS (E)  = 


Using  Eqns.  (3.4)  and  (3.2),  we  obtain 


(3.4) 


00 

i  =  4tt  |  M  |2  f  NX(E)  N2(E+v)  [f  (E)  -  f(E+v)]  dE 

-00 

00 

=  o  /  n1(E)  ng (E+v)  [f  (E)  -  f(E+v)]  dE 

—00 


where  c 


4ir  |  M  I2  Nln(0)  N2n(0)  . 


(3.5) 
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(i)  Normal -Normal  Case  1 


nn 


In  this  case  under  the  usual  conditions  of  v  «  E. 
and  T  «  Ep,  the  current  reduces  to 


‘tin 


cv. 


(3.6) 


Ni  (0)  and  Ng  (0)  are  considered  to  be  constants.  Pig.  (3.1  a) 
illustrates  the  densities  of  states  and  occupation  of  levels  in 
two  normal  metals  near  the  Fermi  surface  at  non-zero  tempera¬ 
tures  and  shows  the  associated  i-v  characteristic. 


(ii)  Normal -Superconductor  Case  ing 

When  one  of  the  two  metals  is  superconducting , 
following  BCS  theory,,  the  density  of  states  then  has  an  energy 
gap  2A^  centered  at  the  Fermi  level.  The  corresponding 
density  of  states  and  i-v  characteristic  are  shown  in  Fig.  (3.1  b) . 
In  this  case 

ins  =  c(v-  -  A^)  ,  I  v  |  > 

(3-7) 

=  0  ,  I  v  |  <  A1 

Also  i  (-v)  =  -  ins(v)  . 
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Figure  3.1 

Energy  diagrams  illustrating  the  density  of  states  near  the 
Fermi  levels  occupation  of  states  and  current-voltage 
characteristic  in 

(a)  normal -normal  case 

(b)  normal -superconductor  case 


(c)  superconductor-superconductor  case. 
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At  T  ^  0°K,  there  is  a  small  current  even  for 
v  <  due  to  the  excited  "electrons"  which  can  tunnel  into 

the  other  metal.  As  the  two  sides  of  the  insulating  layer  now 
look  entirely  different,  we  will  get  a  strongly  temperature 
dependent  current  in  this  voltage  region.  The  energy  gap 
parameter  will,  however,  be  almost  unambiguously  measurable 
despite  thermal  smearing  at  v  %  provided  T  «  A^ . 

(iii)  Superconductor-Superconductor  Case  igs 

The  corresponding  "electron"  density  and  i-v 
characteristic  are  shown  in  Fig,  (3.1  c).  There  are  energy 
gaps  2A1  and  2A^  centered  at  the  Fermi  levels  of  the  two 
superconductors,  Shapiro  et  al  (1962)  have  studied  this  case 
at  different  temperatures. 

At  T  =  0°K,  current  i  will  be  zero  as  long  as 

s  s 

the  applied  voltage  is  less  than  A^  +  A^.  At  v  =  A^  +  A^, 
it  will  increase  sharply  and  will  approach  inn  at  higher 
voltages.  However,  at  a  non-zero  temperature,  some  current 
will  flow  at  the  smallest  applied  voltage  due  to  excited 
"electrons"  above  the  gaps.  This  current  will  increase  until 
v  attains  the  value  Ag  -  ^  (with  Ag  >  &1 ) .  With  further 
increase  in  voltage,  this  current  will  decrease  as  the  same 
number  of  "electrons"  are  available  for  tunneling  as  in  range 


' 


23 


v  <  A2  -  but  into  a  lower  density  of  available  states  in 

the  opposite  superconductor.  This  decrease  in  current  con¬ 
tinues  until  v  is  just  less  than  A1  +  Ag.  As  illustrated  in 
Fig.  (3*1  c),  at  v  =  A^  +  A^,  the  current  rises  sharply 
with  voltage  and  finally  approaches  the  value  i  .  Calcula¬ 
tions  of  Shapiro  et  al  (1962)  show  that  the  singularity  at 
v  =  A2  -  A1  is  actually  logarithmic  and  there  is  a  discon¬ 
tinuous  current  jump  at  v  =  A,  +  A^.  However,  allowance  has 
to  be  made  for  the  factors,  e.g.,  life  time  and  anisotropy 
effects,  neglected  in  the  BCS  theory  and  experimentally  we 
instead  observe  a  cusp-like  peak  and  a  current  jump  spread 
over  a  small  voltage  interval.  However  due  to  piling  up  of 
states  at  the  edges  of  the  energy  gaps  in  both  superconductors, 
the  "k-gT"  smearing  is  of  much  less  importance  in  this  case 
than  in  a  N-S  system. 

B.  Modifications  of  the  Simple  Tunneling  Model 

Harrison  (1961)  has  shown  that  band  structure  in 
normal  metals  cannot  be  observed  through  tunneling  experiments. 
This  is  because  the  density  of  states  is  inversely  proportional 
to  the  group  velocity  of  electrons  in  an  independent  particle 
picture  and  thus  the  density  of  states  in  normal  metals 
completely  drops  out  in  the  tunneling  current  expression 
Eqn.  (3*5)*  In  a  superconductor,  however,  this  independent 


5 


. 
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quasi -particle  model  is  not  true.  In  this  case,  following 
Bardeen  (1961),  the  matrix  element  still  involves  the  group 
velocity  of  the  normal  electrons  and  not  of  the  superconduc¬ 
ting  quasi -particles  whereas  the  density  of  states  involved 
in  the  transition  is  obviously  that  for  the  quasi -particles . 
Since  the  matrix  element  M  is  essentially  unchanged,  the 
only  significant  change  in  the  tunneling  current,  due  to  the 
transition  of  at  least  one  of  the  two  metals  of  the  tunnel 
junction  to  the  superconducting  state,  comes  from  the  density 
of  states  factor  for  the  quasi -particles  in  a  superconductor 
and  it  is  this  alone  that  we  observe  in  the  electron  tunneling 
experiments . 


After  Bardeen Ts  treatment,  Cohen  et  al  (1962) 


described  the  tunneling  of  electrons  through  a  tunnel  junction 
in  terms  of  an  effective  Hamiltonian 


(3.8) 


H  -  H1  +  H2  +  Ht 


where  and  are  exact  Hamiltonians  for  the  metals  1  and 
2  respectively  and  H^  is  an  operator  which  transfers  electrons 
from  one  to  the  other. 


(3.9) 


where  c 


re  c/^and  c,  destroy  and  create  electrons  in  Bloch  states 

ks  k T  s 


of  momentum  1?  and  kT  respectively;  s  denotes  the  spin. 


. 
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Schrieffer,  Scalapino  and  Wilkins  (1963)  treated 
as  a  small  perturbation  on  the  zero  order  Hamiltonian 

Hq  =  Hx  +  H2  .  (3.10() 

Starting  with  zero  temperature,,  if  |  0^>  and 
|  C>2>  are  the  ground  states  of  and  respectively, 

|  0^>  |  C>2>  is  the  ground  state  of  Hq  as  and  are 
assumed  to  commute.  Using  first  order  perturbation  theory, 
the  transition  probability  per  unit  time  for  an  electron  from 
metal  1  to  metal  2  is 


CD. 


1-2 


=  2i t  Z 

a,p 


<0^  |  <P0  |  2  M 


kk's 


y(2)0d) 

kkT  k's  ks 


°i> 


02> 


6(€a+6p-v) 


(3-11) 

where  |  a1>  and  |  P2>  are  the  eigen-states  of  H1  and  H 2 
in  absence  of  perturbation  and  applied  voltage  v.  Ga  ls 
the  energy  difference  between  states  |  a1>  and  |  0-L>  . 
eQ  has  a  similar  meaning. 


Assuming  Mkk,  to  be  constant  in  the  energy  range 
of  interest  (v  «  Ep),  we  can  write  the  transition  rate  and 
hence  the  current  as 

V  (2)  (1) 

/  Nt+  (E)  Nt_  (v-E)dE 
o 


i  oc 


(3.12) 
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where 


and 


-U)  o 

2  I  I  '  Cjj-g  I  01>  |  5(€a-E). 

K] .  CX 


(3.13) 


4^(1)  = 


nI4  (E)  ,  (l  =  1,2)  Is  the  effective  tunneling  density  of 
states  for  adding  electrons  (+)  or  holes  (-). 


Essentially^  electron  interaction  through  phonons 


is  not  instantaneous.  Particle  damping  effects  due  to  phonon 
emission  also  have  to  be  taken  into  account.  In  a  simple  case., 
when  these  factors  are  ignored,  one  obtains  from  Eqn.  (3.13) 
(see  e.g.  Schrieffer  1964) 

Nt+  (E)  =  N(0)  J  E  *  2  1  •  (2-19) 

(E^  -  l\  ) 2 

It  is  to  be  noted  that  NT+  are  just  the  same  den¬ 
sity  of  quasi -particle  states  used  in  the  simple  model  earlier. 

In  real  metals,  Schrieffer  et  al  (1963)  using 
Green’s  function  techniques  have  taken  into  account  a  non- 
instantaneous  potential  for  electron  interaction  through  pho¬ 
nons  and  have  found  a  very  simple  expression  for  NT±  (E) 


(2.20) 


( cf .  Eqn .  (2.19)  ) • 
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A  is  now  a  complex  function  of  energy  E,  the 
imaginary  part  of  which  arises  from  damping  effects  caused 
by  real  phonon  emission.  For  E  »  A,  Eqn.  (2.20)  becomes 

A  2(E)  -  A„2(E) 

Nt±(E)  -  N(0)  [1  +  A - 2—A -  +  .  .  .  ] 

2E  (3.14) 

where 

A(E)  =  A1(E)  +  i  A2(E)  ,  1  =  yTT  . 

This  theory  has  been  further  extended  to  non-zero 
temperatures  by  Schrieffer  (1964). 

We  are  perhaps  now  in  a  position  to  draw  a  clearer 
picture  of  tunneling  phenomenon  in  a  superconductor  which 
will  also  explain  the  complete  absence  of  any  coherence 
effects.  It  is  important  to  note  that  in  the  tunneling  process, 

there  are  two  channels  k  and  Ic1  such  that  kT  is  above  and 

— ► 

k  is  below  the  Fermi  level  with  E^r  =  E^f  .  In  a  N-S  system 
as  shown  in  Fig.  (3.2),  an  electron  in  a  state  Jet  below  the 
Fermi  level  in  the  normal  metal  can  tunnel  into  a  state  TcTt 
above  the  Fermi  level  in  the  superconductor  with  a  probability 
uf ,  which  is  the  probability  that  the  pair  state  (Jc’t  -JeTl) 
is  unoccupied  (see  Eqns .  (2.12)  and  (2.13)).  It  can  also 
tunnel  into  a  state  kt  where  states  JtTT  and  kT  are  related 
by  *  The  Probability  of  tunneling  through  this 
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Figure  3.2 

Two  possible  final  states  for  a  given  initially  occupied 
state  let  in  a  N-S  system. 
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2 

channel  is  u^  which  is  the  probability  that  the  pair  state 

— — >■ 

kt  ,  -kj,)  is  unoccupied.  Since  ekT  =  -e—  ,  u|  equals 
2 

vkt  •  Thus  the  total  probability  that  the  electron  can  enter 

either  state  is  u^,  +  u|  =  u^T  +  v^T  =  1.  This  is  also 

true  when  holes  are  to  be  injected  and  we  then  have 
2  2  2  2 

vkr  +  vk  ~  vk!  +  uk?'  =  ^or  total  probability  of 
injecting  a  hole.  This  is  the  physical  origin  of  cancellation 
of  coherence  effects  and  tunneling  current  may  be  calculated 
by  summing  over  all  states  in  the  superconductor  above  (or 
below  )  the  Fermi  level  only. 

Further  when  an  electron  leaves  the  normal  metal, 
a  hole  is  created  in  it  giving  an  excitation  energy 
e  =1  e,  |  for  this  metal.  When  this  electron  occupies  the 
state  ltrT  in  the  superconductor,  an  excitation  energy 

=  E^,  is  given  to  the  superconductor.  This  can  occur  only 
if  v  =  |  |  +  Ek?  .  Clearly  the  current  will  begin  at  v  =  A. 

In  a  S-S  system,  shown  in  Fig.  (3-3) *  an  electron 
from  It  T  in  superconductor  1  can  tunnel  into  lo T  T  or  It"  T 
and  the  total  probability  that  it  can  enter  either  state  is 
unity  as  before.  An  electron  kt  with  energy  =  -€r-  could 
also  tunnel  to  the  same  final  states.  The  total  probability 

that  an  electron  is  available  for  tunneling  in  any  of  these 

— >  —> 

initial  states  is  that  pair  states  (Tc  T  ,  -Tcf  )  and  (kt  ,  -k  4-  ) 


< 
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Figure  3-3 

Two  possible  final  states  for  each  of  the  two  initially 
occupied  states  let  and  kt  in  a  S-S  system. 
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are  occupied,,  i.e.  v^.  +  =  1.  Thus  in 

both  initial  and  final  states.,  coherence  effects  are  not 

present .  At  T  =  0  K,  as  in  a  N-S  system.,  tunneling  current 

is  possible  only  when  v  ^  v  .  where  v  .  =  A  +  A^. 

x  mm  mm  1  2 

cl  i 

The  dynamic  conductance  (^)  of  a  tunnel  junction 

is  obtained  by  taking  the  derivative  of  Eqn.  (3.12).  The 
d  i 

ratio  of  (-—)  ,  when  one  of  the  metals  in  a  tunnel  junction 

®  d  i 

is  superconducting,  to  (-^r)  when  both  metals  are  normal  is 

n 

called  the  normalised  dynamic  conductance  a  of  the  tunnel 
junction.  For  a  N-S  system  at  T  -  0°K,  one  obtains  from  Eqn. 
(3.12)  using  Eqn.  (2.20) 

(di/dv)s  Nt(F(v) 

(di/dv)n  =  Nt(D(0) 

Comparing  this  expression  with  Eqn.  (2.20),  one  notes  that 
measurement  of  a  will  provide  detailed  information  of  the 
superconducting  tunneling  density  of  states  that  in  turn  is 
a  function  of  the  gap  parameter  which  is  energy  dependent. 

C .  Double  Particle  Tunneling  in  Superconductors 


v 


(v2  -  A2(v)f 


(3.15) 


Taylor  and  Burstein  (1963)  observed  currents  between 
two  superconductors  at  low  temperature  in  excess  of  the  usual 
current  associated  with  single  particle  tunneling  (referred  to 


. 
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as  spt  hereafter)  discussed  above.  These  excess  currents 
were  found  to  be  independent  of  temperature  and  polarity  and 
appeared  with  a  sharp  jump  at  v  =  A1(t)  and  at  v  =  A2(t) 

An  explanation  has  been  proposed  by  Schrieffer  and  Wilkins 
(1963)  in  terms  of  a  double  particle  tunneling  (henceforth 
called  dpt)  mechanism  involving  a  pair  dissociation  or 
recombination . 


In  a  dpt  process.,  two  particles  which  form  a  con¬ 
densed  pair  in  either  initial  or  final  state  are  transferred 
across  the  insulating  layer.  In  a  typical  process,  shown  in 
Fig .  (3.4  a),  which  begins  at  v  =  A^(T),  two  electrons  are 

extracted  from  superconductor  After  tunneling  to 

(2) 

they  recombine  adding  a  condensed  pair  to  S'  .  Thus  no 

(2) 

excitations  are  created  in  S'  A  In  sum  total,  this  leaves 
behind  two  quasi -particles  in  and  the  minimum  energy 

required  for  the  onset  of  this  process  is 


2v 


E 


^ )  +  e/1)  =  2An (T) 


k 


1 


k, 


1 


(3.16) 


Similarly  a  second  process,  shown  in  Fig.  (3-4  b), 

involves  removal  of  a  condensed  pair  from  which  tunnels 

(2) 


to  form  two  quasi -particle  states  and  in  S 

(1) 


excitations  are  created  in  S 


Thus  no 
The  threshold  energy  for 


this  mechanism  is 
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Figure  3**4 

'  r 

Two  mechanisms  contributing  to  double  particle  tunneling 
processes  between  two  superconductors. 


la) 


0 

0 


lb) 


Unoccupied  State 
Occupied  State 
Pair  State 


34 


2v  =  E^)  +  E^)  =  2A2(T)  .  (3-17) 

This  explanation  leads  to  temperature  insensitive,  polarity 
independent  excess  currents  as  observed. 

These  processes  involve  a  second  order  matrix 
element  and  should  result  in  far  smaller  currents  than  those 
observed  for  spt  processes.  However,  as  Schrieffer  and 
Wilkins  (1963)  point  out,  the  effective  thickness  of  the 
insulating  layer  for  the  dpt  and  spt  processes  could  be 
different  leading  to  significant  dpt  currents.  Such  an 
abnormality  in  the  apparent  thickness  of  the  insulating  layers 
involved  in  dpt  processes  has  been  supported  by  experiments 
due  to  Adkins  (1963) » 

D.  Phonon  Density  of  States 

It  was  pointed  out  in  Chapter  I  that  the  phonon 
density  of  states  is  reflected  as  structure  in  the  dependence 
of  ct  on  v.  At  energy 

v  =  Aq  +  0^  ( 3 . 1 8 ) 

where  Aq  is  the  gap  parameter  at  the  edge  of  the  energy  gap 
and  ^  is  the  transverse  or  longitudinal  phonon  peak  energy, 
excited  electron  states  may  readily  decay  by  phonon  emission 


... 


' 
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to  the  increased  density  of  available  states  at  the  gap  edge. 

As  pointed  by  Schrieffer  et  al  (1963)*  this  strong  phonon 
emission  rate  produces  a  sharp  decrease  in  the  cr:v  curve  due 
to  a  rapid  increase  of  A^  and  slight  decrease  in  (see 
Eqn.  (3.14)).  Thus  a  sharp  drop  of  cr  indicates  a  peak  in 
the  density  of  states  of  phonons. 

E .  Recombination  Processes  in  Lead 

Scalapino,  Wada  and  Swihart  (1965)  and  Swihart, 

Scalapino  and  Wada  (1965)  have  recently  calculated  the  energy 

dependence  of  the  gap  parameter  at  non-zero  temperatures  for 

Pb  which  is  a  strongly  coupled  superconductor.  They  have  used 

the  same  phonon  distribution  and  Coulomb  pseudopotential  as 

used  by  Schrieffer  et  al  (1963)  which  is  in  good  agreement  with 

tunneling  experiments  of  Rowell  et  al  (1963) .  However,  in 

Pb,  Scalapino  et  al  (1965)  estimated  the  electron  phonon 

2 

coupling  constant  a.  for  longitudinal  phonons  to  be  double 

2 

the  coupling  constant  a^_  for  transverse  phonons,  to  be 
compared  with  a  constant  value  of  coupling  constant  employed 
by  Schrieffer  et  al  (1963). 

Scalapino  et  al  (1965)  have  calculated  that  there 

may  be  an  anomaly  in  the  effective  tunneling  density  of  states 

of  Pb  at  T  %  T  .  This  is  explained  in  terms  of  removal  of 

c 
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an  electron  from  the  normal  metal  and  a  quasi -parti cle  from 
the  superconducting  Pb  and  creation  of  a  ground  state  pair  in 
Pb  with  emission  of  a  phonon.  At  non-zero  temperatures,  there 
are  a  relatively  large  number  of  excited  quasi -particles  and 
high  density  of  states  at  the  gap  edge  in  Pb.  An  electron 
can  be  injected  from  a  normal  metal  into  Pb  across  an  insula¬ 
ting  layer  and  can  then  combine  with  an  excited  quasi -particle 
and  form  a  ground  state  pair.  A  phonon  is  emitted  in  this 
process.  Because  of  the  high  density  of  states  at  the  gap 
edge  and  the  requirements  of  energy  conservation ,  this  pro¬ 
cess  is  most  probable  when 

V  +  Ao(T)  =  (3.19) 

where  oj  is  a  peak  in  the  phonon  density  of  states  ,  A  (T) 
is  the  gap  parameter  at  the  gap  edge  at  temperature  T.and 
v  is  the  energy  of  the  injected  electron.  Notice  the  differ¬ 
ence  between  the  mechanisms  involved  in  Eqns .  (3.18)  and  (3.19). 

According  to  Scalapino  et  al  (1965)*  this  recombina¬ 
tion  process  causes  a  new  negative  peak  in  the  A2(E,T),  as 
shown  in  Fig.  (3.5)*  which  is  eventually  reflected  in  the 
effective  tunneling  density  of  states.  This  additional  struc¬ 
ture,  due  to  recombination,  is,  however,  subject  to  thermal  smear¬ 
ing.  They  estimate  that  the  maximum  observable  structure  should 
occur  near  T/T  =0.95  around  an  energy 


v  =  3.9  mv . 
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Figure  3 • 5 

Computed  real  and  imaginary  parts  of  the  complex  gap  function 

A  =  A^  +  iAg  as  a  function  of  energy.  The  dashed  curves 

are  for  T  =  0*  and  the  solid  curves  for  a  reduced  temperature 

T/T  =  0.98  (after  Scalapino  et  al  (1965)). 
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CHAPTER  IV 


EXPERIMENTAL  METHOD 


A.  Sample  Preparation 

The  tunnel  junction  samples  were  prepared  in  the 
manner  described  by  Adler  (1963)  and  by  Rogers  (1964)  except 
for  variation  in  the  oxidation  process.  Salient  features  of 
the  sample  preparation  are  reported  here.  The  distinct  stages 
of  preparation  of  a  metal/metal  oxide/metal  sandwich  are 
shown  in  Fig.  (4.1). 

(i )  Base  Layer 

99 . 999^  pure  aluminum  was  used  as  the  metal  for  the 
base  layer  in  this  work.  A  cleaned  glass  slide  was  placed 
close  behind  an  appropriate  mask  in  the  evaporator.  The 
evaporations  were  made  by  heating  in  vacuum  a  0.060"  diameter 

tungsten  wire  carrying  an  aluminum  charge.  All  evaporations 

-5 

were  started  at  a  pressure  of  5  x  10  mm  of  mercury  or  less. 
The  base  layer  was  approximately  1  mm  wide  and  100  to  500  A° 
thick. 

Aluminum  has  been  most  commonly  used  as  the  base 
layer  metal  because  of  the  ease  with  which  it  may  be  evapora¬ 
ted  and  oxidized  to  yield  a  suitable  barrier  layer.  It  also 
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Figure  4.1 

Various  stages  of  sample  preparation: 

A.  A  base  layer  A  was  vacuum  deposited  onto  a  cleaned  glass 
slide . 

B.  A  barrier  layer  A  0  was  formed  by  oxidation. 

x  y 

C.  A  cover  layer  B  was  vacuum  deposited  across  the  barrier 
layer  to  form  a  cross  strip. 

D.  Leads  were  soldered  at  the  ends  of  the  films. 


A 
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becomes  superconducting  at  low  temperatures  and  thus  provides 
a  possible  S-S  system.  Other  metals  that  may  be  used  for 
base  layer  preparation  are  mainly  Sn^  Mg  and  Pb  (cf.  Giaever 
et  al  1962)  all  of  which  become  superconducting  at  low 
enough  temperatures. 

Magnesium,,  on  the  other  hand,  can  provide  a  normal 
base  layer  in  a  tunnel  junction  (cf.  Rogers  1964).  One  of 
the  difficulties  in  evaporation  of  magnesium  was  that  evapora¬ 
ted  magnesium  covered  the  whole  glass  slide  even  when  the 
mask  was  held  only  2  mm  away  from  the  slide  (cf .  Adler  1963) • 
Films  of  magnesium  with  well  defined  edges  have  been  evapora¬ 
ted  on  the  glass  slides  during  the  course  of  this  work  by 
sudden  heating  of  magnesium  ribbons  approximately 
1.5”  x  0.125"  x  0.006"  in  size.  A  cover  slide  was  mounted  on 
top  of  the  experimental  slide  to  prevent  magnesium  atoms 
reflected  from  various  portions  of  the  evaporator  from  reaching 
the  experimental  slide.  However,  magnesium  films  of  any 
desired  thickness  could  still  not  be  prepared. 

(ii )  Barrier  Layer 

The  base  layer  of  aluminum  was  oxidized  in  an  oven 
at  nearly  200°C.  In  one  specimen  reported  here,  preheated  air 
was  blown  over  the  glass  slide  while  it  was  in  the  oven.  The 


. 
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temperature  of  the  oven  was  maintained  at  nearly  200°C  in 
this  case  also.  The  water  vapour  content  of  the  air  near 
the  specimen  is  possibly  higher  in  the  latter  case  and  this 
seemed  to  accelerate  the  oxidation.  Oxidation  time  in  all 
cases  was  about  one  minute. 

(iii)  Cover  Layer 

A  strip  of  high  purity  lead  was  vacuum  deposited 
across  the  barrier  layer.  The  evaporation  was  done  in  vacuum 
by  heating  a  molybdenum  strip  carrying  a  charge  of  lead. 

(iv)  Lead  Connections 

Finally  a  wire  was  attached  to  each  of  the  four 
ends  of  these  two  strips  using  pure  indium  solder.  Fine  copper 
wires  coated  with  50:50  lead-tin  solder  were  used  for  the 
leads.  The  coating  on  these  wires  becomes  superconducting 
at  low  temperatures  and  this  reduces  the  heat  input  to  the 
specimen.  This  method  of  connecting  leads  after  evaporating 
the  cover  layer  allows  the  oxidation  to  be  carried  out  above 
the  melting  point  of  indium. 
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(v)  General  Considerations  in  Sample  Preparation 

Resistance  and  capacitance  of  the  tunnel  junction 
were  measured  with  a  R-C  bridge  described  by  Rogers  (1964) . 

The  special  feature  of  this  bridge  is  that  it  cannot  subject 
the  tunnel  junction  to  more  than  50  mv.  Four  terminal  measure¬ 
ments  of  the  resistance  of  the  tunnel  junction  were  also  made. 
The  specimens  used  in  the  present  work  had  resistances  in  the 
range  from  100  to  500  0  (ohms)  measured  between  the  metallic 
films  and  from  this  the  barrier  layers  were  estimated  to  be 
a  few  tens  of  A°  thick.  The  resistance  of  the  evaporated 
films  was  always  less  than  10  ft. 

The  thickness  of  the  barrier  layers  reported  by 
Giaever  and  Megerle  (1961)  and  by  Rogers  (1964)  depends  on  a 
number  of  variables.  The  oxidation  rate  strongly  depends  on 
the  oxidation  temperature  and  the  atmospheric  humidity  or  the 
residual  H^O  vapour  in  the  system.  Elevated  temperatures  and 
an  increase  in  humidity  facilitate  oxidation  as  was  well 
borne  out  by  the  present  work.  The  oxidation  rate  in  a 
confined  atmosphere  of  oxygen  at  100  to  400°  C  was  found  to  be 
lower  than  in  open  air  at  the  same  temperature.  The  effective 
barrier  thickness  is  also  dependent  on  the  nature  of  the  cover 
layer,,  evaporation  temperature  of  the  cover  layer  and  finally 
the  rate  of  evaporation  of  the  base  layer  itself.  We  have 
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found  during  the  course  of  the  present  work  that  a  tunnel 
junction  with  almost  any  value  of  resistance  could  be  prepared 
by  suitably  controlling  these  parameters. 

The  characteristics  of  the  specimens  prepared  at 
room  temperature  generally  changed  with  time.  Stable  specimens, 
whose  resistance  and  capacitance  remained  almost  constant  at 
room  temperature,  were  prepared  by  oxidizing  the  base  layer 
in  an  oven.  The  elevated  temperatures  in  the  oven  perhaps 
help  in  achieving  equilibrium  surface  characteristics  of  the 
base  and  the  barrier  layers  which  gave  stability  to  the  specimen. 

For  small  voltages  v,  the  current  i  through  the 
barrier  layer  was  found  to  be  proportional  to  the  applied 
voltage  across  it.  The  current  increased  approximately  expo¬ 
nentially  with  voltage  at  higher  voltages  which  is  consistent 
with  the  results  obtained  by  Fisher  and  Giaever  (1961) . 

Figs.  4.2(a)  and  (b)  show  the  i:v  characteristic  of  a  typical 
tunnel  junction. 

The  barrier  layer  does  not  have  a  uniform  thickness. 
This  is  supported  by  results  of  Fisher  and  Giaever  (1961) 
and  of  this  laboratory  (cf.  Rogers  1964).  Measurable  strengths 
of  tunneling  currents  due  to  dpt  processes  have  also  been 
explained  in  terms  of  non-uniform  thickness  of  the  barrier  layer 
(cf.  Chapter  III).  A  few  fine  metal  bridges  piercing  the 
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Figure  4.2 

Current -voltage  characteristics  of  a  typical  tunnel  junction 
at  (a)  low  voltages  (after  Fisher  and  Giaever  1961), 

(b)  higher  voltages  (after  Adler  1963)- 
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barrier  layer  have  also  been  reported  by  Giaever  (i960)  and 
by  other  workers.  This  results  in  a  slight  current  at  vol- 
tages  well  within  the  gap  even  when  one  or  bothythe  metallic 
layers  of  the  junction  are  superconducting,  A  tunnel  junction 
having  a  number  of  such  bridges  will  have  many  discontinuities 
in  its  i:v  characteristic  as  observed  by  Rogers  (1964). 

The  selection  of  a  good  tunnel  junction  was  rather 
difficult.  As  soon  as  it  was  prepared*  its  resistance  and 
capacitance  were  measured.  Only  those  stable  specimens  which 
had  resistance  lying  in  between  100^2  and  500H  were  chosen. 

A  very  thin  barrier  layer  corresponding  to  a  resistance  less 
than  100£2  may  rupture  at  low  temperatures  .  The  higher  limit 
of  the  specimen  resistance  was  partly  set  by  our  measuring 
equipment  and  by  the  fact  that  the  resistance  of  the  specimen 
increased  at  lower  temperatures.  The  primary  characteristics 
of  a  good  specimen  are  the  presence  of  non-linear  i:v  charac¬ 
teristic*  little  increase  of  resistance  with  decrease  in 
temperature  and  sharp  transitions  from  a  N-N  system  to  a  N-S 
or  a  S-S  system  consistent  with  the  theory.  These  character¬ 
istics  confirm  the  predominance  of  tunneling  current  over  any 
other  conduction  mechanism  and  the  absence  of  metallic 
bridges.  The  latter  will  particularly  obscure  any  current 
due  to  dpt  processes.  Unfortunately*  the  existence  of  these 
characteristics  can  be  completely  confirmed  only  when  the 
specimen  has  been  cooled  to  liquid  helium  temperatures. 
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Figure  4.3 


Schematic^  showing  the  method  of  specimen 
specimen  chamber.  The  drawing  is  nearly 


mounting  and  the 
twice  the  actual  size. 
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B.  Specimen  Mounting 

After  a  good  specimen  was  obtained,  it  was  quickly 
mounted  in  a  small  gas  tight  chamber  illustrated  in  Fig.  (4.3)* 
After  the  specimen  was  mounted,  the  chamber  was  evacuated 
through  a  pumping  tube  and  was  then  filled  to  one  atmosphere 
with  helium  gas  at  room  temperature. 

C .  Cryostat  Description 

The  Al-Al  0  -Pb  specimens  were  studied  at  7°K  and 
x  y 

at  approximately  0.3°K.  Hence  a  He^  cryostat  was  necessary 
for  this  work.  At  the  lower  temperature  both  the  metals  of 
the  tunnel  junction  were  superconducting  and  their  energy 
gaps  would  be  fully  developed. 

The  cryostat  has  been  described  in  detail  by 
Adler  (1963)-  Fig.  (4.4)  illustrates  the  experimental  chamber. 
0.3°K  is  roughly  the  lowest  attainable  temperature  in  this 
cryostat.  Electrical  connections  to  the  specimen  chamber 
were  made  with  lead-coated  copper  wires  0.0025"  in  diameter. 

The  copper  wire  and  helium  in  the  specimen  chamber  served  as 
the  thermal  link  between  the  specimen  and  the  He^  chamber.  A 
germanium  crystal  was  attached  to  the  He^  chamber  to  record 
the  temperature  of  the  specimen. 
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Figure  4„4 


s~j-  d-’uCK AA/' 

Schematic^ showing  the  experimental  chamber  of  the  cryostat 
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(i )  Operating  Procedure  of  the  Cryostat 

In  an  actual  run,  the  outer  can  was  soldered  to 
the  plate  M  after  connecting  the  specimen  chamber  to  the  He 
chamber  C  (cf.  Pig.  4*4).  This  joint  at  M  was  tested  for 
leaks  both  at  room  temperature  and  at  liquid  air  temperature. 
Standard  procedure  for  cooling  to  liquid  helium  temperatures 
was  used  except  for  the  following  special  precaution. 

It  was  noted  that  the  pressure  reading  obtained  in 
the  room  temperature  part  of  the  pumping  line  to  the  outer 
can  is  a  very  poor  indication  of  the  amount  of  the  exchange 
gas  remaining  in  the  experimental  chamber.  Large  amounts  of 
helium  gas  still  remained  in  the  experimental  chamber  even 

i 

-6 

when  the  pressure  in  the  pumping  line  was  5  x  10  mm  of  Hg. 
This  gas  could  not  be  pumped  out  of  the  experimental  chamber 
during  a  period  of  several  hours  when  the  latter  was  at  liquid 
helium  temperatures.  It  provides  a  thermal  link  between  the 
He^  chamber  and  the  outer  can  and  temperatures  lower  than 
about  0.6°K  could  not  be  obtained  under  such  conditions.  The 
helium  leak  detector  was  connected  to  the  vacuum  system  for  the 
experimental  chamber  to  serve  as  a  monitor  for  the  amount  of 
helium  gas  present  in  the  experimental  chamber.  The  tempera¬ 
ture  inside  the  experimental  chamber  was  then  allowed  to  rise 
to  nearly  20°K  with  very  little  liquid  helium  in  the  outer 
dewar.  Thereafter  the  leak  detector  reading  fell  by  a  factor 


. 

. 
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3 

of  nearly  10  in  a  period  of  nearly  six  hours.  The  outer 
dewar  was  filled  with  liquid  helium  only  when  the  exchange  gas 
was  thus  completely  pumped  out.  Liquid  helium  temperatures 
could  be  reached  without  admitting  exchange  gas  to  the  experi¬ 
mental  chamber  at  all  below  liquid  air  temperatures.  However, 
it  was  found  that  more  time  was  required  to  reach  liquid 
helium  temperatures  by  this  method. 

The  temperatures  between  0.3°K  and  7°  K  were  ob- 

q 

tained  with  the  help  of  the  heater  mounted  on  the  He^  chamber 

q 

and  partly  by  regulating  the  pumping  speed  in  the  He  system. 

D.  Temperature  Calibration 

The  temperatures  reported  in  this  work  are  those  of 
q 

the  He  chamber  which  was  in  thermal  equilibrium  with  the 

q 

specimen.  The  temperature  of  the  He  chamber  was  inferred 

q 

from  resistance  measurements  of  a  Ge  crystal  mounted  on  the  He 

chamber.  Four  terminal  resistance  measurements  of  this  crystal 

were  achieved  by  using  an  isolating  potential  comparator  of 

the  type  described  by  Dauphinee  (1953)*  It  is  essential  to 

obtain  a  calibration  curve  for  this  thermometer .  Wires  of 

high  purity  lead,  aluminum  and  cadmium  were  wound  on  a  thin 

q 

copper  cylinder  which  was  mounted  on  the  He  chamber.  The 
resistance  of  the  Ge  thermometer  was  measured  at  the  transition 
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temperatures  of  these  metals.  The  transitions  from  the  normal 

to  the  superconducting  state  for  these  metals  were  determined 

by  using  a  four  terminal  resistance  measurement  with  a  circuit 

described  by  Rogers  (1962).  A  fourth  calibration  point  was 

obtained  by  admitting  exchange  gas  to  the  experimental  chamber 

4 

and  opening  the  needle  valve  to  connect  the  He  chamber  to 
helium  gas  at  atmospheric  pressure.  The  vapour  pressure  of 
the  boiling  liquid  helium  was  noted  from  the  mercury  manometer. 

E .  Measurement  of  i:v  Characteristics 

The  i:v  characteristics  of  the  tunnel  junctions  were 
measured  with  a  curve  tracer  described  by  Rogers,  Adler  and 
Woods  (1964)  which  permits  sensitive  four  terminal  resistance 


calculated  at 


measurements 


different  voltages  v  across  the  specimen  from  these  measure¬ 
ments.  Calculations  were  carried  out  with  IBM  7040  computer. 
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CHAPTER  V 

EXPERIMENTAL  RESULTS  AND  DISCUSSION 

A.  Temperature  Measurements 

Two  lead  specimens  with  A1  as  the  base  layer  were 
studied.  Accurate  measurement  of  the  temperature  was  essential 
for  this  work  as  we  were  looking  for  the  additional  peak  in 
the  density  of  states  in  lead  at  temperatures  close  to  its 
transition  temperature  T  .  The  transition  temperatures  of  Pb, 
viz,  7*193 °K  (cf.  Franck  and  Martin,  1961)*  Al,  viz,  1.196°K 
(cf.  Cochran  and  Mapother,  1958)  and  Cd,  viz,  0.52°K  (Martin,  1961) 
were  used  to  calibrate  a  germanium  resistance  thermometer. 

Another  calibration  point  was  obtained  at  the  boiling  point 
of  the  liquid  helium. 


The  above-mentioned  calibration  points  were  success¬ 
ively  fitted  to  the  following  interpolation  equations  relating 
the  thermometer  resistance  and  its  temperature.  These  equations 
are  respectively  due  to  Clement-Quinnell  (1952)  and  Clement 
(cf.  White  1957). 

T"1  -  aUnR)"1  +  b  +  c  InR  (5-1 ) 


and 

T  =  A  log10R/(log10R  -  B)2  .  (5.2) 


■ 
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The  resistance  temperature  curves  obtained  through  these 
equations  are  shown  in  the  Pig.  (5*l).  A  reversal  of  the  curve 
obtained  from  Eqn.  (5«l)  at  higher  temperatures  is  similar  to 
that  reported  by  Rogers  (1964)  for  a  220  ohms  ’speer1  carbon 
resistor.  It  indicates  that  Eqn.  (5.1)  is  not  suited  particu¬ 
larly  to  these  thermometer  materials  over  this  range  of  temperature. 

Prom  these  four  calibration  points.,  three  (or  two) 
were  used  in  finding  the  unknowns  in  Eqn.  (5.1)  (or  Eqn.  (5.2)) 
and  the  remaining  one  was  used  to  check  the  accuracy  of  the 
calculated  temperature  at  the  corresponding  resistance  of  the 
germanium  thermometer.  We  also  had  access  to  the  resistance- 

■3f 

temperature  curve  obtained  by  Dr.  Franck  for  a  similar  ger¬ 
manium  thermometer  which  was  calibrated  by  vapour  pressure 
measurements  over  boiling  He^  and  He4  between  3.88°K  and  0.6l  K. 
Using  this  calibration  curve  the  error  could  be  determined  in 
any  temperature  calculated  with  one  of  the  interpolation  equations. 
The  temperatures  calculated  from  Eqns .  (5.1)  and  ,(5*2)  were 
found  to  be  much  different  from  the  correct  temperatures. 

A  quadratic  equation 

log10R  =  a  +  ^  +  -5-  (5.3) 


It  is  a  pleasure  to  acknowledge  my  gratitude  to  Dr.  J.P.  Franck 
for  providing  me  this  calibration  curve  for  the  germanium 
thermometer. 
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Figure  5«1 

Resistance  temperature  calibration  curve  for  the  germanium 
thermometer  using 

A.  Clement  -  Quinnell  Equation  . 


B.  Clement  Equation  . 


•  Calibration  point* 
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was  tried  and  found  to  be  much  more  successful.  The  resistance- 

temperature  curve  obtained  with  it  is  shown  in  Fig.  (5.2). 

The  uncertainty  in  a  temperature  calculated  through  this 

equation  is  estimated  to  be  ±  0.05°K  in  the  range  from  7.2°K  to 

2 . 5°K  and  was  about  Yf>  at  lower  temperatures.  The  error  in  a 

temperature  close  to  any  one  of  the  calibration  points  is 

much  less  than  this  maximum  uncertainty  and  since  most  of  the 

observations  in  the  present  work  were  taken  close  to  T  for 

c 

lead.,  the  accuracy  of  this  calibration  was  certainly  adequate. 
The  change  in  thermometer  characteristics  after  thermal  cycling 
between  room  and  liquid  helium  temperatures  noted  by  Rogers 
(1964)  in  case  of  the  ’speer’  resistor  was  observed  not  to 
occur  with  the  germanium  element,  but  the  calibration  was 
checked  at  the  liquid  helium  boiling  point  after  each  set  of 
measurements  to  guard  against  changes.  The  current  used  in 
the  measurement  of  resistance  of  the  germanium  thermometer  was 
2  |ra  above  1°K  and  0.5  |ia  below  1°K. 

B.  Recombination  Effects  on  the  Tunneling  Density  of  States 

in  Lead 

Lead,  being  a  strongly  coupled  superconductor  with 
a  relatively  high  transition  temperature,  is  one  of  the  most 
studied  elements  in  tunneling  work.  Giaever  (i960  a)  and 
Giaever  et  al  (1962)  were  first  to  study  the  variation  of 
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Figure  5*2 

Resistance  temperature  calibration  curve  for  the  germanium 
thermometer  using  Eqn.  (5*3). 
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normalized  dynamic  conductance  with  voltage  across  the  specimen 
as  well  as  the  phonon  effects  in  lead.  Rowell  et  al  (1963) 
reported  detailed  experimental  results  for  lead  along  with  the 
theoretical  explanations  by  Schrieffer  et  al  (1963)  and  Scalapino 
and  Anderson  (1964).  Scalapino  et  al  (1965)  have  extended  these 
calculations  to  non-zero  temperatures.  According  to  their 
calculations,  recombination  processes  make  themselves  apparent 
in  the  effective  tunneling  density  of  states  at  temperatures 

close  to  T  for  lead. 

c 

Fig.  (5.3)  shows  the  enlarged  traces  of  two  Av:v 
photographs,  one  taken  when  the  cover  layer  metal  was  supercon¬ 
ducting  and  the  other  when  it  was  in  the  normal  state.  Av  is 
a  measure  of  non-linearity  of  the  specimen  at  voltage  v  (cf. 
Rogers  et  al  1964) .  These  traces  have  been  projected  onto 
graphs  prepared  to  fit  the  actual  oscilloscope  calibration. 

One  of  the  fundamental  things  of  interest  to  us  is 
the  value  of  the  energy  gap  in  lead.  Fig.  (5.4) shows  the  trace 
of  the  actual  i-v  characteristic  for  specimen  Pb-46  at  0.34°K. 
Following  Giaever  et  al  (1962),  the  energy  gap  value  is  deter- 

k 

mined  by  extrapolating  the  linear  part  of  the  jump  in  the  i-v 
curve  of  the  specimen  to  the  voltage  axis.  As  explained  on 
pages  22  and  23,  AB  in  Fig . (5*4)represents  2(A^  +  A^)  where 
A^  and  A^  are  the  energy  gaps  in  the  two  superconductors  forming 
the  specimen.  The  negative  resistance  region  in  the  i-v 
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Figure  5.3 


Enlarged  graphs  from  photographs  of  Av:v  characteristics  for 

(a)  a  N-N  system 

(b)  a  S-S  system 

over  the  same  voltage  range  across  the  specimen. 
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Figure  5  A 

Graph  of  i-v  characteristic  at  the  energy  gap  for  Pb-46 
projected  from  the  corresponding  photograph  at  T  =  0.34°K. 
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characteristic  of  the  specimen  is  not  apparent  in  this  curve 
due  to  the  current  scale  chosen  here.  Rogers  (1964)  found  the 
value  of  2A(0)  for  A1  to  be  (0.42  ±  0.02)  mv  in  this  laboratory 
as  compared  to  0.40  mv  due  to  Giaever  et  al  (1962).  As  our 
method  of  specimen  preparation  was  quite  similar  to  that  of 
Rogers ,  we  shall  use  the  value  of  2A(0)  for  Al  due  to  Rogers. 

This  yields  the  energy  gap  value  2A  =  2.68  mv  for  Pb.  As  this 
work  was  carried  out  at  0.34°K  which  is  much  less  than  the 
transition  temperatures  of  both  Pb  and  Al*  the  correction  to 
get  the  energy  gap  value  for  Pb  at  absolute  zero  is  negligible. 

The  energy  gap  value  2A(0)  for  Pb  is  thus  estimated  to  be 
(2.68  ±  0.06)  mv  from  the  present  work.  This  is  in  good 
agreement  with  the  energy  gap  values  for  Pb  reported  in  the 
literature  (cf.  Rogers  1964).  Giaever  and  Megerle  (1961) 
had  found  the  energy  gap  in  Pb  to  be  2Ap^  =  (2.68  ±  0.06)  mv 
at  1°K. 

In  order  to  carry  out  investigations  of  the  normalized 
dynamic  conductance  vs,  voltage  for  the  Pb  specimen  at  tempera¬ 
tures  close  to  the  transition  temperature  Tc  of  the  thin  cover 
film  of  Pb,  it  was  necessary  to  determine  Tc  accurately.  Fig. (5-5) 
illustrates  some  of  the  i:v  presentations  for  Pb-47  specimen  at 

different  temperatures  in  the  neighbourhood  of  T  .  We  have 

c 

already  seen  that  the  i-v  characteristic  of  a  N-N  system  is  a 
straight  line  and  that  of  a  N-S  or  a  S-S  system  has  a  flat 
region  at  the  origin  (cf.  Fig.  (3.1))  if  the  energy  gap  in  at 
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Figure  5-5 

Graphs  of  i-v  characteristics  projected  from  the  photographs 
around  zero  voltage  for 

A.  T  =  6.4l°K. 

B.  T  =  6. 23°K. 

C.  T  =  6.09°k. 

D.  T  =  5.78°K. 

Both  i  and  v  are  in  arbitrary  units. 
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least  one  of  the  layers  constituting  the  specimen  is  well 
developed.  At  the  transition  temperature,  this  energy  gap  is 
not  fully  developed  (cf.  Fig.  (2.2))  and  there  is  also  a 
relatively  large  number  of  quasi -particles  excited  above  the 
Fermi  level  in  the  superconducting  layer  of  the  specimen. 
Because  of  these  two  factors ,  only  a  short  flat  region  appears 
at  the  origin  in  the  i-v  characteristic  and  it  grows  as  the 
temperature  is  lowered.  The  appearance  of  the  kink  at  the 
origin  is,  however,  a  definite  indication  of  the  transition 
of  at  least  one  of  the  films  of  the  specimen  from  normal  to 
superconducting  state.  Thus  in  Fig.  (5.5.)*curve  A  shows  a 
N-  N  system  and  curve  B  shows  that  the  Pb  film  has  become 
superconducting  (obviously  there  is  no  question  of  superconduc¬ 
tivity  of  the  A1  film  at  this  temperature).  The  flat  portion 
at  the  origin  gets  larger  at  lower  temperatures  as  shown  in 
curves  C  and  D  in  Fig.  (5*5)  .  From  these  results  we  find  the 
transition  temperature  of  this  Pb  film  to  be  6.32°K. 

It  is  worth  noting  that  results  for  this  specimen 
and  for  all  others  observed  in  this  laboratory  indicate  that 
transition  temperatures  of  thin  films  of  Pb  are  lower  than  that 
of  the  bulk  material.  In  contrast  to  this,  the  transition  tem¬ 
peratures  of  thin  films  of  A1  and  In  are  higher  than  those  of 
the  corresponding  bulk  material.  This  is  an  interesting  result 
and  it  will  be  the  subject  of  further  investigations  in  this 
laboratory . 
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Figs.  (5.6  a,b)  and  (5. 7  a)  illustrate  the  cr:v 
characteristics  for  the  Pb  specimen  Pb-47  at  three  different 
temperatures,  /'the  c r:v  curve  at  a  reduced  temperature 
T/T  =  0.95  as  calculated  by  Scalapino  et  al  (1965)  is  shown 
in  Fig.  (5.7  b)  for  comparison.  The  agreement  between  the 
theoretical  and  the  experimentally .observed  curves  is  remark¬ 
able.  As  expected  from  their  calculations ,  structure  due  to 
recombination  effects  is  most  evident  in  the  energy  region 
near  E  =  3-9  mv  (cf.  Chapter  III).  Further  the  flat  region 
between  the  two  peaks  in  the  a:v  curves  gets  shorter  as  the 
temperature  decreases.  This  can  also  be  explained  qualitatively. 
As  the  temperature  is  lowered,  there  are  fewer  quasi -parti cles 
excited  above  the  energy  gap.  The  probability  of  recombination 
of  a  tunneled  electron  with  an  excited  quasi -particle  to  emit 
a  phonon  and  form  a  ground  state  pair  will  then  obviously  be 
reduced.  As  the  recombination  rate  decreases,  the  criv  curve 
will  gradually  return  to  its  usual  shape.  No  definite  ex¬ 
planation  for  the  slight  difference  between  the  observed  and 
calculated  values  for  the  maximum  value  of  cr  at  the  gap  edge 
could  be  found.  This  difference  was  present  in  both  the  spe¬ 
cimens  Pb-46  and  Pb-47-  It  may  be  due  to  inhomogeneity  and 
thinness  of  the  Pb  film. 
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Figure  5*6 


Normalized  dynamic  conductance  a  vs.  voltage  v  across  the 
specimen  Pb-47  at  T  ~  T  . 

(a)  T/T0  =  0.89  . 

(b)  T/Tc  =  0.96. 


(  a) 


(  b) 
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Figure  5*7 

Normalized  dynamic  conductance  c  vs.  voltage  v  across  the 
specimen  Pb-47  at  T  ~  T  . 

(a)  T/T  fc.O.98. 

0 

(b) A. Calculated  curve  by  Scalapino  et  al  (1965)  at  T/Tc  =  0.95  . 

B.  Observed  results  for  Pb-46  at  T/T  =  0.97* 
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Figure  5.8 

Normalized  dynamic  conductance  a  vs.  voltage  v  across  the 
specimen  Pb-46  showing  phonon  structure  at  0.34°K. 
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C .  Phonon  Density  of  States  for  Lead 

The  observations  for  ct:v  for  specimen  Pb-46  were 
also  taken  at  0.34°K  and  the  results ,  which  are  shown  in 
Pig.  (5.8),  are  in  excellent  agreement  with  the  published 
work  of  Rowell  et  al  (1963)  and  the  results  reported  in 
Rogers1  thesis  (1964).  The  origin  on  the  voltage  axis  has 
been  chosen  at  (A^p  +  Ap^)  so  that  the  phonon  structure  energies 
can  be  read  directly  from  the  voltage  axis.  The  BCS  curve 
for  A(E)  =  constant  =  Aq  is  also  plotted  for  comparison. 

The  structure  in  the  experimental  curve  at  energies  corresponding 
to  (v  -  -  Ap^)  equal  to  4.5  mv  and  8.5  mv  is  well  resolved 

and  corresponds  to  the  predominant  phonon  energies  for  lead. 
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CHAPTER  VI 


CONCLUSIONS 


A.  Results  for  Lead 

A  peak  in  the  effective  tunneling  density  of  states 
in  Pb  due  to  recombination  effects  has  been  found.  It  is 
reflected  as  an  additional  peak  near  3*9  mv  in  the  normalized 
dynamic  conductance  curve  for  a  tunnel  junction  containing 
lead.  Our  results  are  in  agreement  with  the  calculated  results 
of  Scalapino  et  al  (1965).  The  variation  of  the  recombination 
effect  with  temperature  has  also  been  studied.  The  recombina¬ 
tion  effect  has  been  found  to  be  most  prominent  just  below  T  . 
The  slight  difference  between  our  experimental  values  and  the 
calculated  values  of  Scalapino  et  al  (1965)  for  a  at  the  gap 
edge  could  not  be  decisively  explained. 

B.  Thermometry  Results 

An  interpolation  formula  with  three  unknowns  has 
been  found  suitable  for  a  germanium  thermometer  for  temperature 
calibration  in  the  range  of  temperatures  from  nearly  7>2°K  to 
0.52°K.  Temperatures  calculated  with  this  equation  are  accurate 
to  ±0.05°K  above  2.5°K  and  to  nearly  1%  below  2.5°K. 
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C .  Suggestions  for  Further  Work 

During  the  course  of  the  present  work  the  following 
experiments,  which  could  provide  more  definite  information, 
were  thought  to  hold  promise  for  future  work: 

csjj.  ttU' 

(i)  A  study  of  the/transition  temperature  with  the  thickness 

of  Pb  film.  This  should  reveal  information  on  the  mechanism 
involved  in  the  change  of  transition  temperatures  of  the 
films  from  that  of  the  bulk  value.  We  have  found  that  the 
transition  temperature  of  Pb,  a  strongly  coupled  super¬ 
conductor,  in  thin  films  is  lower  than  its  transition  temp¬ 
erature  in  bulk.  In  comparison,  the  transition  temperature 
of  Al,  a  weakly  coupled  superconductor,  is  higher  in  thin 
films  than  its  transition  temperature  in  bulk.  Thus  the 
difference  in  the  transition  temperature  of  the  film  and 
of  the  bulk  material  may  depend  on  the  strength  of  the 
electron-phonon  interaction. 

(ii)  A  search  for  recombination  effects  and  a  study  of  the 

TKis 

phonon  density  of  states  in  Hg.  is  interesting  because 

it  is  a  strongly  coupled  superconductor  like  Pb.  It  has 
been  only  slightly  investigated  due  probably  to  difficulties 
in  preparation  of  a  Hg  specimen. 
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(iii)  Use  of  Mg  as  a  base  layer  metal.  Mg  will  remain  as  a 
normal  metal  at  all  lower  temperatures.  The  electron 
density  of  states  in  a  normal  metal  does  not  appear  in 
the  a:v  curve  (cf.  Chapter  III)  and  hence  a  N-S  system 
with  a  base  layer  of  Mg  will  be  helpful  in  the  study  of 
the  phonon  density  of  states  of  the  superconducting 
member  of  the  tunnel  junction. 
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